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Comment on "Intermittent Synchroniza- 
tion in a Pair of Coupled Chaotic Pen- 
dula" 

In pI, a number of supposedly novel and surprising 
features were observed in a system composed of two pe- 
riodically driven and asymmetrically coupled pendula. In 
particular it was claimed that 'permanent synchroniza- 
tion ... does not occur except as a numerical artifact'. It 
was suggested that this might be related to the particular 
type of coupling. In this comment I want to point out 
that some of these claims cannot be maintained. The 
synchronization in this system is precisely of standard 
blow-out type [HJ^. The observed intermittency is ex- 
actly the on-off intermittency well known from the syn- 
chronization of multifractal chaotic attractors y-|§l . 

The system studied in [1] is described by 0m + 6m/Q + 
sin 6*^ = rcos(ni), 9^ + ds/Q + s'mds = T cos{ilt) + 
c[sin 9s — sin 9m] ■ The subscripts m and s stand for mas- 
ter and slave. When they are nearly synchronous, the 
difference 6 = 9m — 9s satisfies the linearized equation 
5 + 5/Q +{l — c)5 cos 6* = 0. My first observation is that 
the same linearized equation would foUw from a sym- 
metric coupling, where master and slave have coupling 
terms ±c/2[sin6'm — sin 0s]. Since the behavior near the 
synchronization threshold is governed by the linearized 
equation, it follows that any eventual abnormal behavior 
cannot result from the asymmetry of the coupling. 

As shown in [EJ , the synchronization threshold is given 
by the condition that the largest Lyapunov exponent 
Ai(c) of the linearized equation is zero. In the present 
case this gives Cc = 0.7948 for the parameter values con- 
sidered in [1] . The eigenvalues of the instantaneous sys- 
tems given in eq.(9) of [1] are irrelevant, except that their 
fluctuations suggest that also the pointwise Lyapunov ex- 
ponents might fluctuate. This is indeed the case. Let us 
consider a finite but large time T and define by Ai and 
A2 the multipliers along the stable resp. unstable man- 
ifold of the linearized equation. Of course they depend 
parametrically on the trajectory 9{t). At the synchro- 
nization threshold, we have (log|Ai|) = rAi(c) — 0, 
where the average is taken over all initial conditions 
9 [to), 9 {Iq)^ 5 {to), 5 [to) with t^ -C 0. Generically we ex- 
pect that (log |A2|)/T = A2(c) < Ai, as is verified numer- 
ically. Therefore we have only one direction in the space 
spanned by [5, 5) along which we must study a possible 
break up of synchronization. 

This break up can occur, even for c > Cc, if Ai fluctu- 
ates and if the system is perturbed by noise pM] . More 
precisely, if this noise is infinitesimal and the fluctuations 
of Ai follow normal central limit behavior, one expects 
intermittent bursts with power behaved distributions of 
amplitudes M and of legths of the locked phase . 

To describe the fluctuations of Ai we use the generat- 
ing function [Q g{z;c) = T"^ log(|Ai|''). The cumulant 
expansion of log|Ai| corresponds to a Taylor expansion 


g{z;c) = z\i{c) + z'^a'^{c)/2 + . . ., where Ta'^{c) is the 
variance of Ai, and contributions of higher order cumu- 
lants are straightforward to compute. The arguments of 
H can now be used straightforwardly to show that am- 
plitudes A = \5\ of the bursts are distributed according 
to P(A) - A^'^-i with g{z = k;c) = 0. Neglecting 
higher order cumulants this gives k = 2Ai(c)/cr^(c). For 
the parameters used in [1], simulations with T — 200 
give o'^(c) = 0.89 at c = Cc, while Ai(c) w 6.1(cc — c). 
The predicted power laws for P(A) are compared to 
numerical simulations in fig.l. In the same Gaussian 
approximation, the distribution for the locking inter- 
vals T is for c < Cc given by the distribution of re- 
turn times to a reflecting wall of a biased 1-d ran- 
dom walk with drift Ai(c) and diffusion constant cr^(c), 
P{t) - r-3/2exp(-rAi(c)V2cr2(c)) §,§. This disagrees 
with the fit in fig. 2 of [1] by the prefactor r~'^/^ which 
gives indeed most of the r-dependence seen in that figure. 
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FIG. 1. Log-log plot of AP(A) for c = 0.795 (threshold), 
c = 0.785 (desynchronized) , and c — 0.805 (synchronized). 
The system was disturbed by noise with level 10~^ . The 
straight lines have the theoretically predicted slopes. 
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